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Abstract. There has been much progress in recent years in understanding 
the existence problem for wave maps with small critical Sobolev norm (in 
particular for two-dimensional wave maps with small energy); a key aspect in 
that theory has been a renormalization procedure (either a geometric Coulomb 
gauge, or a microlocal gauge) which converts the nonlinear term into one closer 
to that of a semilinear wave equation. However, both of these renormalization 
procedures encounter difficulty if the energy of the solution is large. In this 
report we present a different renormalization, based on the harmonic map heat 
flow, which works for large energy wave maps from two dimensions to hyper- 
bolic spaces. We also observe an intriguing estimate of "non-concentration" 
type, which asserts roughly speaking that if the energy of a wave map concen- 
trates at a point, then it becomes asymptotically self-similar. 



1. Introduction 

Let n,m> 1 be integers, and let R^^" be Minkowski space {{t, x) : t £ R, x £ R"} 
with the usual metric gapx^x^ = —dt'^ + dx'^. We parameterize this space by Greek 
indices a,/3 = 0,l,...,n, raised and lowered in the usual manner. 

Let {N, h) be a smooth connected complete m-dimensional Riemannian manifold 
without boundary, and let / C R be an interval. A Schwartz wave map from a slab 
/ X R" in Minkowski space to N is defined to be any smooth map (p : I x R" — + N 
which decays rapidly to some constant (j){oo) Cz N sd, spatial infinity, and whose 
derivatives all decay rapidly to 0.^, and is a (formal) critical point of the Lagrangian 

{d°'(p{t,x),da(f>{t,x))h(^^(^t dtdx. (1) 

The Euler-Lagrange equation for Schwartz wave maps can then be written as 

(0*v)"a„0 = 0, (2) 



1991 Mathematics Subject Classification. 35J10. 

"'^In other words, for each fixed t, <f>(t, x) converges to </>(oo) faster than any negative power of 
and all derivatives of (f> converge to zero faster than any negative power of \x\. We need to 
consider Schwartz wave maps rather than classical wave maps (which are constant outside of in a 
finite union of light cones) because we will need to take (nonlinear) Littlewood-Paley projections 
of these wave maps, obtained via the heat flow for harmonic maps, in our argument. 
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where ((/)*V)" is covariant differentiation on the vector bundle (p*{TN) with respect 
to the pull-back 0*V via (f) of the Levi-Civita connection V on A''. Wave maps are 
of course the Minkowski space analogue of harmonic maps. 

In this paper wc shall restrict primarily to the case when n = 2, m > 2 (the m = 1 
case being reducible to the free (linear) wave equation) and when N has constant 
negative curvature k < 0, thus 

[Vx, Vy]Z - V[x,Y]Z =: R{X, Y)Z - K((y, Z)hX - {X, Z)hY) (3) 

for all smooth vector fields X,Y^Z on A^. Note that k can be normalized to — 1 
by scaling the metric h. A typical example of such a manifold N (with k = —1) 
is the hyperbolic space H™ := 0{m,l)/0{m) = SO{m,l)/ SO{in), where 0{m) 
is the orthogonal group on R™ and 0(m, 1) is the Lorentz group on R "^"j with 
0(to, 1) and SO{m, 1) being the respective orientation-preserving components of 
these groups. More generally, one can quotient by a discrete subgroup of 
0{m, 1) (acting on the left in the usual manner), and as is well known one can 
recover all constant negative curvature m-dimensional manifolds A'' this way, up 
to isomctry and normalization of n. The case A'^ = is of particular interest 
as it arises as a model problem for the Einstein equations under some additional 
symmetry assumptions (see e.g. [1] for a discussion). 

Remark 1.1. With the constant negative curvature assumption, Schwartz maps 
(p become topologically trivial (contractable to a constant map); by compactifying 
infinity, this remark is equivalent to the assertion that any map from a sphere to 
N is topologically trivial. When N is hyperbolic space H™ this claim is clear since 
A^ is contractible; more generally, if A^ is quoticnted out by a discrete group, 
then the claim follows by lifting the map up to hyperbolic space (the discreteness of 
the fibers keeps the monodromy trivial). Later on we will see that these maps can 
in fact be explicitly deformed into a constant map by means of the harmonic map 
heat flow. Note that while A^ itself need not be parallelizable or even orientable, 
this will not particularly concern us since maps into A'' can be lifted up into H"*, 
which is both parallelizable and orientable. Indeed one could reduce to the case 
N = H™ by lifting the wave map up to hyperbolic space, and observing that the 
wave map equation (2) is preserved by this lift. 

In local co-ordinates = on A'', the wave map equation (2) can be written as 

d"d„r = -vU4>)dact^''d''4>'' (4) 

where the F^^ are Christoffel symbols; in particular we see that (j) evolves by a non- 
linear wave equation. It is then natural to consider the Cauchy problem, specifying 
initial data (t>[to\ := {(l){to) , dt(l){to)) at some initial time to G R. Of course, this 
initial data ^[to] must take values in the tangent bundle TN of A'', and decay rapidly 
to (^!>(oo), 0) as X ^ oo; we refer to such data as Schwartz initial data. 

The Cauchy problem for wave maps has been extensively studied (see references); 

we refer the interested reader to the surveys in [23], [30], [50], [59], [69]. For 
arbitrary dimension n and arbitrary targets A^, one can easily show by general 
energy methods (which indeed apply to any nonlinear wave equation, see e.g. [56]) 
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that for any Schwartz^ initial data <j>[to], one has a unique local-in-time Schwartz 
solution to (2) in some spacetime slab / x R", where / is a neighborhood of to, 
attaining the value of ^(oo) at infinity (thus the value of </> at spatial infinity is 
fixed by the wave map evolution). Indeed, energy methods can be used to show 
such smooth solutions can be continued whenever the norm \\Vt,x4'{t)\\H'^/^+^(R'^) 
remains bounded for some £ > 0. In particular, by use of the conserved energy 

E{cl>) = E{m) ■■= I \\'^t,.<i>{t,x)\l dx (5) 

Jar 2 

we see that one has global existence of Schwartz solutions in one dimension; see 
[16], [33], [12], [46], [21]. The one dimensional wave map equation is actually rather 
special, being completely integrable and thus supporting special solutions such as 
homoclinic breather solutions; see [48]. 

Now we turn to higher dimensions n > 2. Since the wave map equation (2) is 
invariant under the scaling 

(6) 

for any A > 0, we see that the two-dimensional case n = 2 is critical with respect to 
the energy (5), as this energy is now scale-invariant (dimcnsionlcss) in this case. For 
higher dimensions n > 2 it is known that singularities can form^ in finite time when 
the manifold N is positively curved, [46] , [53] , [5] , [50] (and there is strong numerical 
evidence to suggest that there is also blowup when n — 2 [18], [2] in this case); when 
n>7 one can also obtain singularities for negatively curved manifolds [5], despite 
such manifolds being well-behaved for other equations (such as the harmonic map 
heat flow). However, one still has global solutions for small data. This was first ob- 
served for initial data sufficiently close to a point [6] or a geodesic [55] , and for data 
which is small in the critical Besov norm ||Vt,a;(?!il|^n/2-i^j|^"^ [66], [67], assuming 

that the manifold N is uniformly smooth. This was achieved by viewing the wave 
map equation as a nonlinear wave equation in (p (as in (4)) and iterating in a certain 
sophisticated function space. These Besov space results appear to be essentially the 
limit of the method; for instance, it is known that such iteration methods cannot 
work if the data is merely assumed to be small in 1| Vt.a;0]| h"/^~'^(R,") ^^^^ [21]i [39], 
[9]). However, this small Sobolev norm global regularity result turned out to be 
achievable upon either a microlocal renormalization of the wave map equation, or 
by passing to a derivative formulation and selecting a special gauge (such as the 
Coulomb gauge). In the case of the sphere = S*™^^, this was achieved for n > 2 
in [64], [65] via the microlocal renormalization method; for the more general class 
of boundedly parallclizable manifolds (which includes the hyperbolic spaces H™), 
this was achieved for n > 5 in [27] (by a hybrid of the microlocal renormalization 
and Coulomb gauge methods) and then for n > 4 in [51], [38] and n > 3 in [31] (by 
use of the Coulomb gauge). For the hyperbolic plane iV = the n = 2 case was 

^Thc smoothness condition has been relaxed substantially, see for instance [26], [29], [21], [66], 
[67] , [68] , [69] and the discussion below; however we shall restrict our attention here to the Schwartz 
category. The rapid decay assumptions can also be removed by finite speed of propagation, though 
possibly at the cost of creating a domain of existence which is not a spacetime slab / X R" . 

•^Despite the presence of these singularities, it is known that global weak solutions exist globally 
in time, see [36], [11], although it is quite possible that uniqueness and energy conservation will 
fail for these solutions. 
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recently treated in [32] (again using the Coulomb gauge) , while for manifolds which 
are uniformly isometrically embeddable in Euclidean space the result was obtained 
for all n > 2 in [68] (together with more precise well-posedness results). While the 
result in [68] does not directly cover the hyperbolic spaces H"*, which cannot be 
uniformly isometrically embedded, it may be possible that the argument can be 
modified to treat this case by first quotienting H™ by a discrete group in order to 
compactify the target, and lifting back to H"* at the end of the argument. 

Thus we now have a fairly satisfactory regularity theory when the critical Sobolev 
norm is small; in particular in two dimensions we have a small energy regularity 
theory for a reasonably large class of target manifolds. Now we turn to the question 
of large energy data at n = 2. In light of the numerical work in [18], [2] we do not 
expect global regularity in the positive curvature case, but it was still conjectured 
(see e.g. [24]) that one has global regularity for large energy data in negative cur- 
vature manifolds such as H"*; note that this would in fact imply an unconditional 
global regularity result for smooth data thanks to finite speed of propagation. In 
the case where the wave map has some U{1) symmetry (either rotation-invariant, 
or rotation equivariant) , such results have already been obtained in [7] (for spheri- 
cally symmetric maps) and [52], [49], [50] (for equivariant maps). The spherically 
symmetric regularity result in fact extends to arbitrary targets [60], [61] but the 
equivariant one seems restricted to sufficiently "negatively curved" targets; see [62] 
(and the numerics in [2], [18]) for further discussion of this issue, which is in partic- 
ular related to the existence of non-constant harmonic maps from S'^ to N, which 
among other things can be used to generate non-trivial finite-energy stationary so- 
lutions to (2)). numerics that strongly suggest blowup occurs for instance when the 
target is a sphere, even for equivariant data. In [62] it is shown that for equivariant 
data, blowup can only occur if the wave map relaxes to a (rescaled) harmonic map 
from to 5^ for a sequence of times converging to the blowup time. 

The conjecture of large energy global regularity for wave maps into H"*, if true, 
would be analogous to other critical large energy global regularity results in the 
literature, notably that of the energy-critical nonlinear wave equation Du = in 
three dimensions (see e.g. [50]), the critical Yang-Mills-Higgs equation in three 
dimensions [19], and the energy-critical nonlinear Schrodinger equation iut + Au = 
\u\'^u in three dimensions (see [14], [4] for the radial case, and [8] for the general 
case). These large energy results are substantially more difiicult than the small 
energy theory, and require at least two new ingredients: 

• A perturbation, local existence, and theory suitable for large energy solu- 
tions, assuming that a suitable spacetime norm (e.g. the L^*^ norm for the 
energy-critical three-dimensional nonlinear Schrodinger equation, or L^L]^ 
norm for the energy-critical three-dimensional nonlinear wave equation or 
Yang-Mills-Higgs equation) is bounded; 

• Some sort of energy non-concentration argument which establishes that 
such a spacetime norm is bounded. 

The energy non-concentration argument generally relies on some sort of mono- 
tonicity formula, often arising from a Morawetz inequality, although in some cases 
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(notably that of the energy-critical Schrodinger equation) such an inequality is not 
sufficient by itself to control the desired spacetimc quantity; one must supplement 
that inequality with additional arguments, such as the induction-on-energy strategy 
pioneered by Bourgain [4] and which also plays a fundamental role in [8]. On the 
other hand, in these types of arguments, the large energy perturbation analysis does 
not require such a Morawetz inequality, and instead proceeds by a variant of the 
standard existence theory, using Strichartz estimates as the primary tool. In order 
to do this, one relies crucially on the fungibility in time of spacetime norms such as 
the ; what we mean by this is that if the Lj^^ norm is large but bounded, then 
one can decompose the time interval into a bounded number of sub-intervals, such 
that the L(*J. norm is small on each interval. This allows one (at least in principle) 
to obtain a perturbation theory for large L^*^ norm solutions from that for small 
norm solutions. This is in contrast with norms such as L'^L^ which are clearly 
bounded by the energy for (say) the energy-critical NLS, but are not fungible in 
time and so cannot be used directly to establish a global regularity result for this 
equation. 

We do not fully resolve either of these two issues here, however we can announce 
some progress on both. More precisely, we are able to establish a large energy 
perturbation theory for wave maps, in the case when the target is a space of constant 
negative curvature, although the actual details of this theory are rather technical 
and will be only partially sketched here. Secondly, we have a simple monotonicity 
formula argument which shows that if the energy of a wave map (with arbitrary 
target) concentrates at a point, then it must be asymptotically self-similar in an 
averaged sense (see Proposition 3.1). 

We now discuss the gauge renormalization (leaving the discussion of the non- 
concentration result to Section 3). For the model problem of constant negative 
curvature, a direct attempt to mimic either the microlocal renormalization proce- 
dure in [64], [65], [27], [68] or the Coulomb gauge renormalization procedure in 
[51], [38], [31], [32] runs into difficulty (although the Coulomb gauge renormaliza- 
tion will work in the abelian case, in which the target is a two-dimensional space 
such as H^, so that the gauge group 50(2) is abelian). More specifically, one be- 
gins to encounter difficulty in large energy in keeping the microlocal gauge change 
approximately unitary, and in the Coulomb gauge one has problems establishing 
uniqueness, regularity, and ellipticity of the gauge in the non-abelian case where the 
target has dimension m > 2 and so the gauge group SO{m) is non-abelian. Also, 
the microlocal gauge has not yet been successfully applied to manifolds such as hy- 
perbolic space, which are not easy to embed inside Euclidean spaces; meanwhile, the 
Coulomb gauge introduces some additional (and somewhat artificial) singularities 
at the spatial frequency origin ^ = (arising from the Hodge decomposition) which 
cause additional technical complications, especially in two dimensions (sec [31], [32] 
for a demonstration of this phenomenon). Finally, while both these gauges arc well 
suited for establishing regularity of a single solution, the task of establishing a per- 
turbation theory - i.e. comparing two nearby solutions - is also not trivial, as one 
must then establish that the gauge transforms themselves have suitable continuity 
properties, which is an extremely delicate issue in the critical setting. 
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In this report wo give a new renormalization procedure, which is based on the 
harmonic map heat flow, which is intrinsic to the geometry and seems a extremely 
natural choice of gauge; we shall dub it the "canonical heat-temporal gauge" or 
"caloric gauge" . This caloric gauge has essentially the same effect as the microlocal 
gauge change (which can be viewed as a discretized, partially linearized, version of 
the harmonic map heat flow renormalization), in that it replaces the nonlinearity 
by a "paraproduct" variant in which at least one derivative falls on a low frequency 
term. However, unlike the microlocal gauge change, the caloric gauge works well 
in large energies and negative curvature targets, essentially thanks to the famous 
result of Eells and Sampson [10] which shows that the harmonic map heat flow for 
those targets is globally smooth and converges to a constant. The caloric gauge 
also is well adapted to perturbation theory, since the harmonic map heat flow also 
enjoys a good perturbation theory. When compared against the Coulomb gauge, 
the caloric gauge has the advantage of renormalizing the nonlinearity into a slightly 
better form (with no singularity near the spatial frequency origin); we shall give a 
heuristic comparison of the two renormalized nonlinearities using Littlewood-Paley 
analysis later on in this report. It is also insensitive to non-abelian behaviour in 
the gauge group, and does not develop difliculties with uniqueness or regularity, as 
long as the target is negatively curved. 

The discussion here will consist of rigorous geometric computations, and non- 
rigorous analytical heuristics; a rigorous treatment of the perturbation theory which 
is informally discussed here will appear elsewhere (and in a much lengthier paper). 

Acknowledgements. The author is indebted to Ben Andrews and Andrew Hassell for 
a crash course in Riemannian geometry and manifold embedding, and in particular 
to Ben Andrews for explaining the harmonic map heat flow. The author also thanks 
Mark Keel for background material on wave maps, and Daniel Tataru for sharing 
some valuable insights on multilinear estimates and function spaces. The author 
thanks Andrew Hassell and the Australian National University for their hospitality 
when much of this work was conducted. The author is a Clay Prize Fellow and is 
supported by a grant from the Packard Foundation. 



2. The derivative formulation 

Here and in the rest of the paper we fix to be a smooth complete Riemannian 
manifold of constant negative curvature k < 0. We now describe the (standard) 
derivative formulation of the wave map equation (2), writing \/t,x(l> in terms of an 
orthonormal frame e{t,x): this formulation was used for instance in [27], [51], [38] 
(or in [17] for harmonic maps). Our discussions here shall be primarily geometric 
rather than analytic. As we shall need this derivative formulation not just for 
the wave map eqiiation, but also for the harmonic map heat fiow (and in fact 
for a concatenation of the two equations), we shall first describe the derivative 
formulation of a general smooth map : f2 — > AT, where fl is an open contractible 
subset of either a Minkowski space R"'^"''"' or a Euclidean space R". In either case 
we will use the co-ordinates to parameterize the domain fl, with the associated 
partial derivative operators da- For this general discussion we will not need the 
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metric structure of fl (i.e. we will not raise and lower indices), although we will 
use the fact that il is flat in the sense that the partial derivatives da commute with 
each other, thus [da, dp] = 0. 

We recall the orthonormal frame bundle Framc(A^) of N , defined as the collection 
of all pairs (0, e) where cj) £ N and e : R™ T^N is an orthogonal transformation 
(or equivalently, an orthonormal basis of T^N). It is geometrically obvious that, 
Prame(-/V) is a smooth principal 0(TO)-bundle of TV, where 0{m) is the orthogonal 
group on R™ (not necessarily orientation preserving). Since is contractible, we 
can lift any map : O — > A'' to a map (0, e) : f2 ^ Pramc(A^) into the orthonormal 
frame bundle. In fact there are multiple such lifts; given any smooth function 
U : ^ 0{m,), we have the gauge transform {4>,e) {(f), ell^^) which transforms 
any lift of to any other lift of indeed, since Prame(A^) is a principal 0{m)- 
bundle, all such lifts can be obtained from a single specified lift by a (unique) gauge 
transform. 

Let ((/), e) -.Vt ^ Frame(A^) be a smooth map. The bundle (j)*{TN) (the puUback of 
the tangent bundle of A'" under (j)) is an m-dimensional vector bundle on Q,, and the 
inverse of the orthonormal frame e can be used to identify this with the trivial 
bundle Q. x R™. In particular, since the derivatives da(j) are sections of (t>*{TN), 
we can define the R'"-valued one-form tpa : U ^ R™ by this procedure, thus 

ipa := e~^da(j) 

or equivalently 

da(p = etpa (7) 



The Levi-Civita connection V on TN induces a pullback connection ^*(V) on 
(j)*{TN), which in turn induces a connection on the bimdlc il x R™ by the above 
identification. This connection is compatible with the gauge group 0(m) (since the 
original Levi-Civita connection is also), and is given by the covariant derivatives 

T>a--=da+Aa. (8) 

where the A„ : O — > o(m) take values in o(m), the skew-symmetric linear transfor- 
mations of R™ (and the Lie algebra of 0(m)) and are defined by 

(</.*V)ae = eAa. (9) 

Prom the Leibnitz rule and (9) we have 

{rV)a{ef) = e(Daf) (10) 

for any smooth scalar function f : CI R"*. Since the Levi-Civita connection V is 
torsion free, we have 

which upon pulling back by (j) becomes 

{<i>*V)ad0(t> = {cl>*V)0cj>daCl>. 

Prom (7), (10) we thus have the fundamental zero torsion identity 

D„V/3 - D^V'a = 0. (11) 
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This identity allows us (at least in principle) to write any of the fields ipa in terms 
of a single field tpp, thus converting the derivative wave map system into something 
resembling a scalar equation. 

Remark 2.1. The tangent bundle TFrame(iV) of Frame(A^) can be canonically iden- 
tified with Frame(A^) x R™ x o{m), since for any (0, e) S Frame(A^), the orthonor- 
mal basis e can be used to parameterize the horizontal space T^N (lifted up to 
T(^.e)Frame{N) via the Levi-Civita connection) as R™, while the vertical space of 
the bundle Frame(A^) at {(p,e) can be parameterized as o(m). The fields {^pa,Aa) 
can thus be thought of as the derivatives da{(l>, e) of the original fields (0, e) with 
respect to this identification. Furthermore, one can show that the Lie algebra 
structure of this R™ x o(TO)-bundle induces an action of the special Lorentz group 
0{m, 1), which then acts transitively on Frame(A''). This in turn can be used to 
verify the assertion made earlier that N is the quotient of H™ by a discrete group. 
We were however unable to exploit this group structure to any particular advantage, 
and shall mostly treat the fields ^, e and ipa, separately. 

The connection has a curvature tensor Faf} : O — > o(m) defined by 

F^p := [D„,D^] = d^A0 - dpA^ + K,^/?]- (12) 

We can use the constant negative curvature hypothesis (3) to compute Faj3. Indeed, 
from (3) we have 

[(</>*v)„, {(b*v)p]z = K{{dfj(b, z)Hdo,4> - (a„<^, z)Hdp4>) 

for all sections Z of (j)*{TN); applying (7), (10), (12) we obtain the curvature 
identity 

Fal3 = Klpa A V/J (13) 

where we use u Av G o(m) to denote the skew-symmetric linear operator on R™ 
defined by^ 

{u A v)w := u{v, w) — v{u, w). 

This identity will allow us (given a suitable fixing of gauge, see below) to recover 
A as some sort of quadratic integral of ifj. 

Remark 2.2. In the Riemann surface case m = 2, the gauge group 0(m) is abelian 
and so the commutator term [Aa, Ag] vanishes in (12). This is in particular useful 
for placing A in the Coulomb gauge. However in the canonical heat-temporal gauge, 
we will not be unduly bothered by this commutator term, indeed in that gauge this 
commutator will have a number of vanishing components and will thus essentially 
be ignored. 

We refer to {4>,e) as the undifferentiated fields and {ilJa,Aa) as the differentiated 
fields. The gauge freedom e eU~^ mentioned earlier affects these fields by the 



The cancellation inherent in this wedge product operation will be largely unexploited, in 
contrast to [32], [31], where this type of "Qij null structure" is needed to compensate for the 
singularity of in the Coulomb gauge. We will, however, rely crucially on the negativity of k 
to ensure that the heat flow converges properly. 
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formulae 

((/.,e)^ ((/.,et/-i); {^^, A^) ^ {U^P^^UA^U-^ + (daUjU-^); Fo,p^UF^pU-\ 

(14) 

In particular we observe that the magnitudes \^a\ and \Fafj\ are independent of 
the choice of gauge. Here we use \Fa/3\ to denote the Hilbert-Schmidt norm, thus 
|F| := tr(Ff 

Remark 2.3. Clearly from (7) and (9) we can obtain the differentiated fields from the 
undifferentiated fields. Now we consider the reverse question, whether one can ob- 
tain the undifferentiated fields from the differentiated fields. To fix boundary condi- 
tions, let us now assume that O is convex, unbounded, and connected at infinity, and 
that e) converges rapidly (in the Schwartz sense, i.e. all derivatives V^^{(f),e) 
decay faster than any polynomial) to some point (^(oo), e(oo)) G Framc(A'') at 
infinity''' This then implies that (V'a,^a) are Schwartz. To reconstruct (</», e) from 
{ipcti Aa) one simply solves the ODE (7), (9) from the given data at infinity. A sim- 
ple calculation using (11), (13), and Frobenius's theorem^ shows that this system 
is indeed solvable; we omit the details. 

To summarize, we now have have two equivalent ways to view a map (0, e) : O — > 
Frame(A'^); either in an undifferentiated formulation, in which (cj), e) can move freely 
(but do not take values in a linear space), or in a differentiated formulation, in 
which the fields (■i/'a, A^) now live in linear spaces but are constrained to obey the 
compatiblity conditions (11), (13). In the Schwartz category these two viewpoints 
are equivalent as long as we fix (0(oo), e((X))) in advance. As we shall be needing 
both the differentiated and undifferentiated fields at various stages our analysis, we 
shall however not take advantage of this redundancy, and instead consider all the 
fields simultaneously. 

2.4. The derivative formulation of the wave map equation. We now special- 
ize to Schwartz wave maps 4> '■ /x R"" N for some compact interval /, which equal 
some fixed point 4>{oo) G N a,t infinity. Given any fixed frame e(oo) in T^^^^^N , 
we can lift (f) to a. map (0, e) : 7 x R" Frame(A/') which converge rapidly to 
((^(oo), e(oo)) at infinity. 

The wave map equation (2) can be rewritten using (7), (10) as 

where wc of coiirse now \isc the Minkowski metric to raise and lower indices. Com- 
bining this equation with the compatibility conditions (11), (13) and the definitions 

^We remark that if 4> decays rapidly to 4){oo) at infinity and e{oo) is any orthonormal frame 
of T^(cx:,)N, then one can easily find a lift (</>, e) : Q —* Prame(Ar) of </> : f2 — > AT which converges 
to (0(oo), e(oo)) at infinity. If N is hyperbolic space H"" this is clear since H*" is parallelizable 
and one can simply pull back an orthonormal frame on H"" which equals e(oo) at <p{oo); more 
generally one can lift the map (j> : Q ^ N to a map <p : Q ^ (retaining the Schwartz property 
at infinity, obtain a frame (<j), e) for the lift, and then push the frame back down to N. 

Alternatively, if S Q, one can use a radial gauge x" {(f>*'V)ae = 0) and propogate (7), (9) 
along radial directions; one can then use (11), (13) and a Gronwall inequality argument to show 
that (V'ai^a) are indeed the derivative fields of {4>,e). Again, we omit the details. 



10 



TERENCE TAO 



(7), (9) we can thus write the wave map equation as the system (cf. [51], [38], [31], 
[32]) 

(0*V)„e = eAo, 
D"Va = (15) 

Note that these equations are both over-determined (the compatibihty conditions 
(11), (13) enforcing a number of constraints between the fields) and under-determined 
(since we still have the gauge freedom (14)). In light of Remark 2.3 it is possible to 
omit the first two equations, and just consider the last three equations as a system 
for the differentiated fields ijja , Aa ; since can then recover the undiff'erentiated fields 
(j), e using Probenius's theorem (and the boundary data (^(oo), e(oo)) as discussed 
earlier; note that we have relied crucially on the constant curvature in order to 
forget about cj) and e. However, we will retain both the differentiated and undiffer- 
entiated fields in our system, as the presence of the undifferentiated fields makes 
the compatibility conditions automatic. 

We refer to the quantity u := D^^/^q as the wave-tension field (as opposed to the 
heat-tension field ~ ^ktpk which will appear later); thus wave maps are precisely 
those maps whose wave-tension field vanishes. 

The system (15) is not yet well-posed in time, even locally, because of the gauge 
freedom (14). There are a number of options available to fix the gauge freedom 
and recover well-posedness; examples of such gauges include the tem,poral gauge 
A(j — 0, the Lorenz gauge da A" = 0, and the Coulomb gauge djAj = 0, where 
Latin indices such as j will run over the spatial indices 1, . . . ,n with the usual 
summation conventions. The latter gauge has been exploited for instance in [38], 
[51], [31], [32]; one advantage of this gauge is that it makes A a quadratic function 
of tp. Indeed, from Hodge theory one can easily verify in the Coulomb gauge that 

A^ = A-\djdjA„) 

= A-\dj{Fj„ + [A^, Aj]) + d^djAj) (16) 

= A-^dj{Kll)j A Va + [Aa, Aj]), 

where A is the spatial Laplacian. Schematically, we thus have A = W^^lijjtjj + AA). 
This identity allows the wave map equation to be quite tractable in high dimensions 
(see e.g. [38], [51] for the n > 4 cases), however the inverse derivative causes 
some difficulties in low dimensions, although they can eventually be surmounted 
(see [31], [32]). Also in the large energy setting the \/~-^{AA) term also becomes 
problematic, indeed it can cause breakdown of uniqueness in the Coulomb gauge 
(though this issue does not arise in the Riemann surface case m ~ 2). Because of 
these difficulties, we shall choose a different gauge, which we shall call the canonical 
heat-temporal gauge or caloric gauge, which will be described later in this section and 
is constructed using the heat flow for the harmonic map equation; it is the covariant 
formulation of the microfocal gauge used in [64], [65], [27], [68], and places A in 
the schematic form A = 7r(V', \/~^tp), where tt is a (nonlinear) paraproduct, defined 
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using the heat flow for the harmonic map equation, which effectively constrains 
the 'V~^ip factor to have higher frequency than the ip factor. This is similar in 
strength to the Coulomb gauge V~^{tpip), but mostly eliminates the "high-high 
interactions" where two high frequency components of ip interact to create a low 
frequency component of A; this is precisely the term for which the operator 
causes difficulties. Also we do not encounter the V~^{AA) term at all. 

As observed in [38], [51], [31], [32], the derivative field tp evolves by a nonlinear 
wave equation. Indeed, from (15) we have 

and 

and hence 

= d"{-A^^p + AM + a^a^jVa 

= 5«(A^Va - - 5^(^"Va)- ' 

Given that A should behave like a quadratic expression of ip, containing an inverse 
derivative (cf. (16)), we thus see that (17) is something like a cubic wave equation 
for 4>- 

As is well known (see e.g. [50]), the wave map equation enjoys a stress energy tensor 

Ta/3 := {da(l),di3(j))h - ]^gaf){d^<i),d^<j))h, 
or in terms of the differentiated fields 

where g is the Minkowski metric; observe that this quantity is indepedent of the 
choice of gauge e. In particular we can define (in standard co-ordinates x'^ = 
t,x^,x^) the energy density Too by 

Too{t,x) ■.= ^\^IJt,^{t,x)\\ (19) 

Because the domain R^"*"^ is invariant under translations in spacetime, Noether's 
theorem yields the conservation law 

5"T„^ = (20) 

which can of course also be verified directly from (15). The conservation of the 
energy 




now follows directly from (20). 

Remark 2.5. In the case of a Euclidean target N = R™ (so k = 0) and the standard 
orthonormal gauge e = dx^, ■ ■ ■ ,dx^, we have tpa = da4>i A^ = 0, and the wave 
map equation (2) then becomes d°'da(j) = 0, while (17) becomes 9"9aV'/3 = 0- Thus 
in this case the wave maps equation collapses to the free (linear) wave equation. 
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2.6. The heat flow equation. We now leave wave maps for the moment and turn 
to the heat flow equation for harmonic maps (or heat flow for short). This equation 
can be written covariantly as 

5,0 = {ry)kdk^ (22) 

where (p{s. x) G is a smooth map defined for some interval [0, S) xR^ for some < 
S < +oo; formally, this is gradient flow for the energy functional ^|Va;0|^ dx. 
Actually, for our applications we will not consider just a single map 0(s, x) evolving 
by the heat flow, but an entire one-parameter family (/)(s, i, x) of such maps, where 
t ranges in a compact time interval /. While the s variable plays the role of time in 
this equation, it is quite distinct from the time variable t appearing in the wave map 
equation, and so we shall refer to s as the heat-tem,poral variable and t as the wave- 
temporal variable to minimize confusion. We use Latin indices j, k to parameterize 
the spatial variable x (and write x to denote both indices, thus for instance 
denotes the vector {Ai, A2)), and Greek indices a, f3 to parameterize the spacetime 
variables t, x, but shall always treat the heat-temporal variable s separately (using 
neither Greek nor Latin indices to parameterize these). In our discussion of the 
heat flow the wave-temporal parameter t plays almost no role, although we will 
need to study the evolution of time derivatives such as dt4> in the s direction. 

In local co-ordinates (j)"", the equation (22) becomes a non- linear heat equation 

= + r^e(<A)5a/9"<^^ (23) 

as such it is clear (e.g. from energy methods) that this equation is locally well 
posed in the Schwartz class, and the flow fixes the value </'(oo) of the map at 
infinity. Indeed, we have the global well-posedness result 

Theorem 2.7 (Global existence of heat flow). If (j) : I x N is Schwartz 

and equals 4>{oo) € N at infinity, then there is a unique Schwartz solution (j) '■ 
R^ X / X R N to (22) with initial data (p{0,t,x) — (p{t,x). Furthermore, (/){s) 
converges to (?i(oo) in the C"^{I x R^; A'') topology as s ^ +00 for any m > 0. 

This result is essentially a special case of the results of Eells-Sampson [10] (see also 
[44], [35], [34], [58]), and follows primarily from the negative curvature of N. We 
will not prove this standard result here, but we mention that a key tool in this 
convergence result is the Bocher identity 

{ds - A)|V^„|2 = k|V^„ a - 2|D,Va|', 

which in conjunction with the diamagnetic inequality \Vx\ipa\\ < iD.iV'al and the 
negative curvature hypothesis k <0 implies that the evolution of \ipa \ is dominated 
by the heat flow: 

(a, -A)|^„| <0. (24) 

We can also write the heat flow equation in differentiated form by lifting (t> to 
{4>, e) : R"*" X 7 X R^ — > Fiam.e{N) as before. We will however limit the gauge 
freedom here by imposing the heat-temporal gauge condition 

(<^*V),e = 0; (25) 
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in other words, we require that e propagates in the heat-temporal direction ds by 
parallel transport. It is clear geometrically (or from the Picard existence theorem) 
that for fixed (j), the frame e{s,t,x) at {s,t,x) is now determined completely by 
the corresponding frame e{0,t,x) at s = 0. We still have the gauge freedom (14), 
but only if the gauge U = U{t,x) is independent of s, i.e. it depends only on 
t, X. We will remove this gauge freedom shortly, however, by selecting a canonical 
heat-temporal gauge. 

Using this gauge e we can introduce the derivative fields ^/^a , V's , , Ag as before, 
although the heat-temporal gauge condition forces = Q (or equivalently = 
ds). From (10) and (7) we see that (22) becomes = ^ki^k- In analogy to (15), 
we can now write the heat flow as the system 

da<t> = eipa 

ds(j) = eips 

{(p*V)se = ^ = 

i^s = DfeVfe (26) 

Dalpp - 'DpTpa = 
dsi>a = T^si>a = 'Dai's 

dsAa = Fsa = Kips A V'a- 

We refer to the quantity ips = Dh^pk as the heat-tension field; it is the direction 

that the heat flow evolution in the s variable deforms </) in (with respect to the 
orthonormal frame e), and should be thought of as a non-linear version of the 
Laplacian applied to the map (p. 

Similarly to (17), the derivative flelds ipa obey a (covariant) heat equation which is 
linear (but with potential and magnetic components). Indeed from (26) and (12) 
we see that 

dstpa = DaVs 

= Dc.DfcV'fc , ^ 

(27) 

= ^kTialph + Fak->Pk 

= DkDktjjc, + Kilpa A Vfe)'^i'fc; 

this equation can be viewed as a linearization of the heat flow equation around (p. 
Observe that the heat-tension field tps also obeys the same equation: 

dsips = 'DkIik^l-'s + n{ips Aipk)i'k- (28) 

As a first approximation, one can thus think of ipa{s,t) as being a non- linear ana- 
logue of e^^ipa{0, t), where e*'^ is the propagator for the free (linear) heat equation. 
Similarly for the heat-tension field ijjs- As is well known, heat operators can be 
used as a substitute for Littlewood-Paley operators; thus one can think of ipa{s,t) 
heuristically as a nonlinear restriction of ipaiO,t) to frequencies |^| < s~^^^. The 
heat-tension field tps is similar, but the presence of the derivative in the formula 
%ps = Dkipk suggests that this quantity is more localized to the annulus |^| ~ 
than the baU |^| < 5-^2. 
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Remark 2.8. The assumption k < will ensure that the equations (27), (28) will 
not blow up as s — > +oo; this is basically because of (24). 

As mentioned above, wc still have a gauge freedom to rotate e by an arbitrary 
gauge U{t,x) independent of s. However, we can remove this freedom by fixing a 
boundary condition at s = +oo. 

Theorem 2.9 (Existence of canonical heat-temporal gauge). Suppose : 7x ^ 
A'' is Schwartz and equals (/>(oo) G N at infinity, and let </> : R"*" x / x be 
the heat flow extension given by Theorem 2.7. Let e(oo) be an orthonormal frame 
in T0(oo)iV. Then there is a unique lift (?i, e) : R"*" x 7 x R^ — > Frame(A'') which 
is sm,ooth, obeys the heat-temporal condition (25), and which converges uniformly 
to (^(oo), e(oo)) as either s — > +00 or a; —> 00. Furthermore, Aa^s) and ipa{s) 
converge uniformly to zero as s ^ +00. 

We will not prove this theorem in detail here, but follows from decay estimates on 
■0 as s ^ +00 which ultimately stem from (24), which ensure that one can specify 
boundary data at s = +00 in a well-posed manner. 

We refer to the gauge e constructed in this manner as the canonical heat-temporal 
gauge or caloric gauge associated to (j) (and to e(oo)). Geometrically, it is con- 
structed by dragging the constant frame (0(oo), e(oo)) back via parallel transport 
from s = +00 to s = by reversing the heat flow; it can also be viewed as the 
unique solution to (26) with boundary conditions 



Note that we are specifying data at both s = and s = +00. As we shall explain 
shortly, this gauge is the analogue of the microlocal gauge in [64], [65], [27], [68] in 
the differentiated setting, and will be used here in place of the Coulomb gauge that 
was used in [38], [51], [31], [32], for reasons which will also be discussed shortly. 

Note that we can now compute A from ip via the fundamental theorem of calculus 

and (26): 



and in particular we can compute the connection Aa{t,x) on the boundary s = 
by the formula 



(0, e) = ((/)(oo), e(oo)) when s 
{ipcnAa) = when s = +00. 



-|-oo 



(29) 




(30) 




(31) 
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If we continue the heuristic that t]j{s,t,x) is a non-hnear version X), 
then the above expression can then be viewed as a non-hnear analogue of a para- 
product 7r('(/'a;(0), V~^V'a:(0)), whcrc the paraproduct tt is a "low-high" paraproduct 
that restricts the frequency of the V^r^Va;!^) factor to be larger or comparable to 
that of the Va;(0) factor. To see this informally, we restrict s to the dyadic re- 
gion 2~^^ < s < 2"^^+-'^, and then heuristically view ipx{s,t) as P^yApx{0,t), a 

Littlewood-Paley type projection to frequencies < 2^^, while xpsi-s.,t)T)k'>Pk{s,t) is 
roughly like 2^Pjj,^j;(0, i), a projection to frequencies ~ 2^. The expression (31) 
is then heuristically like J2]<iiP<'k''^xiQ,t)){2^^P^^^x{0)), whence the claim. This 
paraproduct A = 7r(V'x(0), V~^Vx(0)) is favorable as the inverse derivative is 
always on the highest frequency factor; this should be compared with the Coulomb 
gauge, which gives a formula for A of the form A = ^^^i^x'^Px) (see (16)), which is 
of similar strength to the paraproduct in (31) except for high- high interactions of ip, 
for which the inverse derivative is in an unfavorable location; this is not particularly 
harmful in high dimensions (as the arguments in [38], [51] show) but begins to cause 
significant problems in two and three dimensions. In [31], [32] the difficulties arising 
from this unfavorable placement of the inverse derivative are overcome by using the 
overdetermined nature of the wave map system to extract additional cancellations 
from the expression '^x^i''Px'>Px), and by quite delicate multilinear estimates. We 
will face similar difficulties, but they will be milder because of the more favorable 
location of the inverse derivative. 

Remark 2.10. We now perform some heuristic manipulations that connect this 
gauge to the microlocal renormalization used for the sphere C R*" in [64], 

[65] (and can also be used to connect to the similar procedures used in [27], [68]). 
Here we think of the wave map <j) as a column vector with m entries. We choose an 
orthonormal frame e = (ei, . . . , e^) taking values not on the tangent space of the 
sphere, but rather in the ambient Euclidean space R™. In the standard orthonormal 
frame (which we call e"), the analogues of ipa and A^ are 

Ipa ■■= da(l); Aa := da4) A (J). 

In a more general frame e = Ue^, where U takes values in the rotation group 0{m), 
we have 

V'a := U-^da(j); A„ := U-\da<p A 4>)U - U'^dJI. 
The heat-temporal gauge condition = thus becomes 

dsU = (9,0(s) A <P{s))U{s) 
where dsip can be defined using the heat flow equation 

= A0 + (;6|V(/.|^ 

While this equation can blow up for large energies, for small energies it converges to 
a constant. In the region s ~ 2"^^^, we can then heuristically treat (j){s) as the low 
frequency projection P^j^(/)(0), and ds(t> ~ as the medium frequency projection 
s~^P]j^(0). Thus the equation for U is then approximately 

sd,u^[p^m^p<km]u{s). 

Comparing this with the scheme used in [64] (see also [65], [27], [68]) we see that 
the canonical heat-temporal gauge construction is nothing more than a continuous 
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version of the discrete microlocal gauge used in those papers, with the Littlewood- 
Paley operators in the ambient EucHdean space being replaced by the more intrinsic 
operators generated by the harmonic map heat flow. Note that in the heat flow 
formulation, the matrices U are automatically orthogonal, whereas in the discrete 
proccdiu-e in [64], [65], [27], [68] one only obtains approximate orthogonality, and 
only when the energy is assumed small. Indeed for large energies it seems unlikely 
that one can replicate this procedure, due to topological obstructions (Schwartz 
maps from to 5™"^ are not necessarily topologically trivial). 

Remark 2.11. The equation (30) can be used to recover A from ipa. One can 
similarly use (26), (30) and the canonical heat gauge condition ipa{oo,t,x) = to 
recover ipa from the heat-tension field ijjs'. 



Thus one can (at least in principle) recover ipa from tps by an iteration process, 
and then by (30) recover Aa from ijjg and ijja. Thus one can view (26) as a scalar 
non-linear heat equation (28) for the heat-tension field tpg We will thus often seek 
to use the above identities to write the other components of ■0 and A in terms 
of this "dynamic variable" ijjg in order to make the system (34) resemble a scalar 
equation (and thus make more amenable to solution by iterative methods); this 
corresponds to the "dynamic separation" strategy used in [31], [32]. Indeed one can 
view the heat-tension field tps as a sort of^ continuous Littlewood-Paley resolution 
of (/), which is the scalar variable in the undifferentiated wave map equation (2). 
One consequence in particula-r from (32) that the functions ijjct are (modulo a cubic 
lower order term '^a) the gradient of a scalar function ^ i'sis) d-s. This fact 
will be crucial for us, as it means that all of the highest order non-linear expressions 
involving the tpa contain null form structure (the lower order terms will turn out 
to be manageable via Strichartz estimates and will not require null structure). 
In particular we will seek to recover null forms such as the trilinear expression 
^1 5" 02 9a 03, preferably with 0i having "higher frequency" than 02 or 03, in the 
spirit of [65]: in our licjat-flow formulation the s parameter will be a proxy for this 
concept of frequency®. A similar strategy was carried out in [31], [32], starting not 
from the heat flow but instead from the equation D„^^ — D/j^a = to conclude 

'''a more precise heuristic is that if s ~ 2^^'^, then i/is --^ ~P]i4'- 

®More precisely, the quantity corresponds to the Fourier notion of frequency magnitude, 

as can be seen by inspecting the Fourier symbol of the heat operator e'"^ ; thus larger values of s 
correspond to lower frequencies and vice versa. 




(32) 



where 'J' is the cubic correction term 
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that "0 has small curl and licncc (by Hodge theory) is approximately a gradient. 
However the use of Hodge theory once again introduces inverse derivatives V""'^ 
into the equation, which makes the non-linearities more difficult to control than 
would otherwise be necessary, and so we shall rely on (32) instead of Hodge theory 
to write ij} as an approximate gradient. 

2.12. Concatenating the wave map and heat flow equations. From The- 
orems 2.7 and 2.9 we see that the system (26) can be solved, and placed in the 
canonical heat- temporal gauge (29), given any Schwartz map ^ : 7 x — > A'' (as 
well as an orthonormal frame e(oo) in T^(^co)^- In particular we can solve this 
equation for any Schwartz solution to the wave map equation (15). This gives rise 
to the fields (j), e, Aa, ipa, ^s, tps on R"*" x 7 x R^ solving the system of equations 
and boundary conditions 



ds4> 


= eijJa 

= eips 




(<^*V)„e 


= eAa 




{(f>*V)se = A, 


= 






= when s = 





V's 


= ^ktpk 






= 




dsi^a = 'Dsi'a 


= DaVs 




Fafj 


= Kljja A Ipp 




(^sAct — Fsa 


= Kljjs A Ipa 






— {(f){oo), e((X) 


)) when s = +oo 




= when s = 


-|-oo. 



(34) 



Of course, many of these equations are redundant, and this system is highly overde- 

termined, nevertheless it will be convenient to retain all of these equations as we 
will need the full structure of this system at various points in the argument. 

Remark 2.13. If the target N were Euclidean space N = R™ with the standard 
frame, then we have V'a = Vs = ^4>j A^ — Ag — 0, and the heat flow equation 
(22) then becomes ds4> = At/), while (27) becomes dgtpa = ^ipa- Thus in this 
case the heat flow collapses to the free (linear) heat equation. It is clear that (34) 
similarly collapses to an evolution which is the wave equation in the x, t variables 
and the heat equation in the x, s variables, with the two flows commuting with each 
other. For non-Euclidean targets, the same statement is true to top order; thus all 
fields evolve in the s direction by a nonlinear heat equation and in the t direction 
by a nonlinear wave equation. 

Remark 2.14. One can use the transitive group action of 0{l,m) of Frame(iV) to 
place 0(oo) and e(oo) wherever one pleases. However there will be no need for us 
to perform such a normalization. 

Remark 2.15. We caution the reader that the heat-temporal variable s scales like 
twice the dimension of space, in contrast to the wave-temporal variable t, which 
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scales like one dimension of space. Indeed, the scale invariance of (34) is given by 
((A, e){s,t,x) 



1ps{s,t,X 
Fal3{s,t,X 



(</',e)(^,|,f) 



1 



s t X , 



-(V.„,A„)(-,-,-) 



(35) 



If (- - -) 

A3 '"^A2'A' A^- 



Note that the wave map equation D^ipa = only holds at the boundary s = of 
the slab R"*" x J x R^. Unfortunately, this equation docs not commute with the 
heat flow equation (even in the completely integrable one-dimensional case n = 1), 
so we do not obtain this equation in the interior of the slab. However, we can assert 
that the wave map equation holds "approximately" in the interior region. Indeed, 
if we define the wave-tension field u hy u := D^^/i^, then we have (from various 
equations in (34) and covariant versions of the Leibnitz rule) 

= K(^,A^„)^" + D"D„Va 

= K(Dk^k A Va)V'" + D"D„Dfe^/>fe 

= KiBki'k A Va)^'" + D"(F„feVfe) + D"DfeD„^fc 

= ^[(DfcVfe A Va)V'" + D"((V'a A Vfc)V'fe)] + D"DfeDfeVa 

= K[(DkA A V'a)^'" + (D^i^a A i^k)A + i^a A T>"i^k)A + i^a A (35) 

+ F„feDfe^° + Dfc(F„fcV") + DfcDfeD«^„ 
= K[(Dfe^fc A Va)^" + {uA '4)k)'(l^k + ii^a A T)k^l^")^IJk + (V'a A Vfe)D"V'fe 

+ 2(^„ A ^fc)DfcV'° + (Dki^a A Vfe)^i'" + (^a A Dfe^fc)V"] + DfeDfeU 
= DfcDfcM + «;(zt A 'ipk)^k + A DfcV^")V^fe + 3(V'a A ipk)'Dk->p" + {Dk->pa A ^fc)^"] 

= DfcDfcU + k{u a V'fe)V'fc + 4K(V'a A 4>k)'Dk'4>"- 

Thus the wave-tension field it propogates in the s direction by a covariant, inho- 
mogeneous heat equation. Also, from (34) we have w = on the boundary s = 0. 
Thus we expect u to be somewhat "small" for s > as well. 

Just as the wave-tension field obeys a heat equation in s, the heat-tension field 
ips = D/;V'fe obeys a wave equation in t: 

a«5„V« = 9"(D„V.-^aV'.) 

= a«a«^„-a«(^„V.) 

= a«(D«-A«)^„-a«(^,V.) 

= a««-a«(A«Va)-s"(A„^.). 
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One can use (36) and (37) to design a (local-in-time) iteration scheme for solving 
(34) in the large energy setting. The rigourous details of this scheme will appear 
elsewhere, but the informal ideas are as follows. We first use (36) propagate the 
wave-tension field u in the s direction (using values of 'ilJt,x and At^x obtained from 
previous iterates), using the boundary condition u = at ,s = 0. Then one uses 
(37) to propagate the heat-tension field tps forward in the t direction (using the 
value of u just obtained, and values of Aa and tps obtained from previous iterates), 
using the initial data of at the initial time t = to. Finally, one uses (32) and (30) 
to compute i/'t.x and At^^ from tpg (and using values of V't,x obtained from previous 
iterates), using the boundary condition Vt.x = ^t,x = at s = -|-oo. This scheme 
may appear convoluted, but it is necessary to do this in the large energy setting 
in order that the iterative algorithms to compute each of the fields ipt,xi ^t,x) 
eventually involve some integration in time (so that we can take advantage of time 
localization). In the small energy case there are simpler schemes available^, but the 
convergence of such schemes in the large energy case becomes problematic because 
of the possibility of an infinite number of iterations backwards and forwards in the 
s variable (which cannot be made to converge by localizing t). 

Remark 2.16. An overly simplified model of the system (34) can be obtained by 
assuming the approximation ipa « for some function $ (cf. (32), (11)), and 
assuming that the nonlinear heat flow behaves like the linear heat equation (and 
in particular behaves somewhat like a family of Littlcwood-Paley operators); then 
the equation for $ on the boundary s = is roughly of the form 

□$« o((Pk,*)(5aPk.^)(^"^k3*)); 

ki,k2,k3:ki>min(k2,k3) 

thus $ evolves by a nonlinear wave equation with non-linearity of the form $3q,$9"$, 

but with at least one of the derivatives falling on a low frequency term. This is 
essentially the type of renormalized equation obtained for the wave map flow when 
the target is a sphere [65] or an isometrically embedded manifold in Euclidean 
space [68] , and it is a good heuristic to keep in mind for the local theory (including 
the large data perturbation theory). For comparison, the corresponding heuristic 
equation for the Coulomb gauge formulation of the wave map equation would look 
something like 

« V-MVx$(5a$)]a"cE. 

(see [51], [38], [31] and especially [32]), which is a similar type of equation but has 
an additional divergence arising from "high-high" interactions inside the term 
that must be controUed^^. 

^For instance, one can propogate iptjx forward in t on the boundary s = using (17), then 
propagate in s using (27), and then recover A using (30). One can pursue these ideas to obtain a 
new proof of the small energy regularity results, at least in the constant negative curvature case, 
which is close in spirit to the arguments in [65] and particularly in [68] , but we will not do so here. 
One could also abandon the caloric gauge, and use the Coulomb gauge instead; for instance, the 
arguments in [32] can (in principle) be modified to obtain an iteration scheme when the target is 
H^. However the V^^ factors arising from this gauge would require one to prove far more delicate 
multilinear estimates, and in particular it becomes even more difficult to take advantage of the 
time localization than in our current argument. Also the Coulomb gauge has some uniqueness 
problems in the large energy non-abelian (m > 2) case, being in some sense even more non-local 
than the caloric gauge. 

^'^There is an additional null structure on the above trilinear form of "Qij type" which can 
be used to handle this divergence, see [32]. However, in the caloric gauge we neither have the 
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3. Asymptotic local self-similarity 

We now discuss the asymptotic self-similarity of wave maps that concentrate at a 
point. Our arguments here are based entirely on an analysis of the stress energy 
tensor Ta/3, and in particular are independent of the choice of gauge or on the 
curvature properties of the manifold. 

Based in prior experience with large data critical non-linear wave equation regu- 
larity results (see e.g. [50] for some examples), one might hope (in the negative 
curvature case) to establish an energy non-concentration result, for instance show- 
ing that 

lim / Too{t,x) dx = (38) 

i|x|<|t| 

whenever V' is a Schwartz derivative wave map on [—1,0) x ; this for instance 
would be sufficient to establish global regularity for large energy wave maps by 
combining such an energy non-concentration result with a small energy regularity 
result (e.g. the one in [32]) and exploiting a finite speed of propagation argument^^. 

As is well known, one can hope to obtain (38), at least for certain components of 
the energy density Tqo, by contracting the stress energy tensor T^p against some 
well-chosen vector field^^ X^^ in spacetime^^ and then applying Stokes' theorem on 
a truncated backwards light cone such as {|a;| < T < t < — e}; the point is 
that (20) allows us to compute the divergence of TajsX^ as^^ 

a«(T„;3X^) = T^fid'^X^. (39) 

The boundary terms on the cone can be handled by a standard flux decay argument 
and can be considered negligible as a first approximation. The strategy is then to 
choose so that the expression in (39) (as well as the boundary terms TofjX^) 
consists of terms which are either non-negative or small. 

In the case of spherical symmetry or equivariant symmetry, one can demonstrate 
energy decay away from the time axis {x = 0}, not directly by the above strategy, 
but rather by taking advantage of the symmetry assumptions to eliminate the angu- 
lar components of the stress-energy tensor, thus reducing (20) to what is essentially 
a set of transport equations for the remaining components of the stress-energy ten- 
sor in the null directions dt±dr. In particular it is relatively easy (via a Gronwall 

divergence, nor the null structure required to deal with that divergence, and so the analytical 
treatment is somewhat simpler in this gauge. 

^^Actually, there is a slight difficulty in truncating the wave map properly to take advantage 
of finite speed of propagation, as one has to take a little care to ensure that the truncated wave 
map still has small energy and obeys the required compatibility conditions; we will not discuss 
these (minor) technicalities here however. See [68], [69] for some closely related issues. 

^^We shall abuse notation and identify vector fields X = with their corresponding first 
order difTerential operators X = Xl^dp without further comment. 

^"^In related equations such as the semilinear wave equation, one often adds lower order cor- 
rection terms to T^pXl^; however there are no such lower-order terms available for wave maps. 

^^One can of course exploit the symmetry of T^g to replace d"X^ by the more symmetric 
deformation tensor ■k'^I^ := ^{d°'Xl^ + dl^ X"), but we shall not need to do so here. 
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inequality argument, see e.g. [50]) to show in these special cases that 




Mt\<\x\<\t\ 



Too(i, x) dx = for any < A < 1. 



(40) 



If one then applies (39) to vector fields such a,s X = rdr = a; • and applies the 
above estimate, one can easily obtain (time-averaged) decay of the time component 
of the energy: 



Prom this and the negative curvature of N we can then obtain decay of the spatial 
component of the energy also. 

In the absence of symmetry, one cannot hope to apply this type of argument directly. 
First of all, the stress-energy tensor no longer propagates solely in the null directions 
dt±dr, but can now propogate in any timelike or lightlike direction, which seems to 
frustrate any attempt to use Gronwall's inequality to obtain much decay, even very 
close to the light cone. A more serious objection, however, arises from the Lorentz 
invariance of the wave map equation once the symmetry assumptions are removed. 
Indeed, if one could prove the decay estimate (40) away from the time axis for 
arbitrary Schwartz wave maps, then by Lorentz invariance one could also prove a 
similar decay estimate away from any other timelike ray emenating backwards in 
time from the spacctimc origin (0.0). Using two such disjoint timclike rays, it is 
then an easy matter to obtain (38). Thus there is no advantage in excluding the 
time axis in proving (40) for general wave maps; this estimate is as difficult as 
the original estimate (38). In particular, one cannot hope to obtain an estimate 
such as (40) without^^ relying more heavily on the negative curvature of N; this 
in particular rules out the possibility of proving (40) purely from an analysis of 
the stress-energy tensor. Similar considerations hold for (41), basically because the 
vector field dt is not Lorentz-invariant. 

However, a closer inspection of (41) shows that, in light of (40), the content of (41) 
is only new when {t,x) is close to the time axis; thus we only need the relevant 
vector field to point in the direction of dt near the time axis. This now allows for 
the possibility of a Lorentz-invariant estimate, for instance using the scaling vector 
field S := x^da = tdt + x ■ Vx, or more generally S/p'^, where p := ^—x'^Xa = 
(|tp — |a;p)^/^ and fc e R is some parameter. For harmonic maps in three dimensions 
it is well known (see e.g. [43]) that the choice k = 1 (which would give the radial 
vector field dr in three dimensions) yields a useful monotonicity formula of Pohozaev 
type; in light of the algebraic analogy between Euclidean space R^ and Minkowski 
space R^"'"^ it is then natural to try the same thing for wave maps on R^"*"^. This 
indeed works, and gives the following decay estimate: 



^^For instance, the numerical work in [2] strongly suggests that when the target is and 
the initial data is cquivariant, then the energy cein concentrate along the time axis leading to 
blowup. Applying a Lorentz transformation, one can then obtain non-equivariant initial data 
which concentrates on another timelike ray, and in particular (40) fails for this choice of data. 




(41) 
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Proposition 3.1 (Asymptotic local self-similarity). Let tp be a Schwartz derivative 
wave map on [—1, 0) x (not assumed to have any symmetry assumptions). Then 



' ' ,2 ^ dt 



[ I \Us{t.xr dx-^=Q, (42) 

T-,0- |l0g|r|| J\:o\<\t\ t t 



where ips is the component of ip in the scaling direction S, thus 



This estimate should be compared with (41). On one hand, it replaces the vector 
field dt by the variant jS = dt + j-'Vx, which is (conformally) invariant under mild 
Lorentz transformations (i.e. Lorentz transformations with bounded coefficients^^) 
and applies for general wave maps. On the other hand, it requires somewhat more 
averaging in time than is present in (41). Note that conservation of energy only 
allows us to say that the expression in the left-hand side of (42) is bounded, but 
not that it decays to zero. 



We now give a proof of the above proposition. 



Proof [Proof of Proposition 3.1] Let ^ be a smooth derivative wave map on [— 1, 0) x 
R^. We now elaborate the stress energy analysis begun in (39). For any two times 
— 1 < i2 < ii < 0, define the truncated solid cone 

%2,ti] ■■={{t,x):t2<t<ti;\x\ < \t\}. 

This truncated cone has the boundary dKf^ U dK^^^^^^j U dKf^, where dKt is the 
disk 

dKt := {{t,x) : \x\ < \t\} C {t} x R^ 

and is the truncated light cone 

dK^t^^t^) := {{t, x):t2<t<tv,\x\ = \t\}. 

The surfaces dKt can be given Lebesgue measure dx (which co-incides with the 
measure induced from the ambient Minkowski metric). The surfaces dK(t2M) ^re 
null and thus have no canonical measure, however we will give them the artificial 
measure da, defined by 

u da := / u{t, \t\ijj) \t\du}dt 



I 

JdK,t 



(t2.*l) 



for all test functions u, where du is arclength measure on the circle .Si. From 

Stokes' theorem we thus have 

-/ d°'Pa{t,x)dtdx= / Poit2,x)dx+ Pl da- Po{ti,x)dx, 

for any one-form Pa smooth on K[t2,ti]^ where 



PL{t,x) :=Po{t,x)-^Pj{t,x) 



^^We do not expect our estimates to be invariant under extreme Lorentz transformations since 
such transformations tend to increase the energy substantially. 
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is the L := dt — dr component of P. In particular, from (39) wc have 
- [ Tc,f3d"X^ dxdt = [ To/sX^ dx + [ Tl/sX^^ da - j Tq/jX^ 

for any vector field X^ smooth on K\t^^t-i_]- 

We can now derive the standard facts of energy monotonicity and flux decay by 
choosing the vector field X = 9t in (43), obtaining the energy identity 

= / Too da; + / ^lo da - j Too dx. (44) 

Prom (18) one easily verifies that 

on 5-ftr(t2,ti). where y := j^^d^^ — j^^d^^ is the spatial angular derivative, thus 
i^iy — ■pj''/'2 — j^^i' particular, we have energy monotonicity 

/ Too dx< j Too dx for all - 1 < t2 < ii < 0; 

if we then set Eq := ^ Too dx, then we have 



dx 
(43) 



Too dx < Eq for all - 1 < i < 0. (45) 

dKt 

Prom this, (44), and monotone convergence we obtain the bounded flux property 

Tlo da <Eq< oo. 



dK, 



(-1,0) 



In particular we have flux decay 



lim / Tio rf<T = 0. (46) 



Now let A > 4 be a large number, and let — 1/A <t\ <Q. We set t2 := Ati, and let 
r]{t, x) = rj(t) be a time cutoff which equals 1 when ^2/2 < t < 2ti, vanishes when 
t < t2 or t > ti, and is smooth otherwise (so in particular T]'{t) = 0{l/t2) when 
t2 <t < 12/2 and r]'{t) = 0{l/ti) when 2ti <t < ti). We would now like to apply 
(43) to the vector field 

p 

where p := ^—x°'Xa = — |a;p, as the left-hand side of (43) will then give 

something very similar to the expression in (42). Unfortunately this vector field 
is singular on the light cone and so (43) does not directly apply. To resolve this 
problem we shall mollify the above vector field slightly. Let < £ < 1 be a small 
parameter, and define the modified co-ordinates x°' and p by 



5p := x° + eti; x^:=x^; p = \J —x'^Xa 
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for j — 1,2; thus the x" co-ordinates are just the x" co-ordinates but with the 
origin shifted upwards to (— eii, 0). In particular p is now smooth on all of K^t^^f^]. 
We now apply (43) with the vector field 

Using the identities 

P 

where g is the Minkowski metric, we see that 

gax/3 ^ ^{x"xf^ - g'^^x'x^) + pO'^d^TJ. 

Contracting this against T^^ using (18) (and the fact that g^^Qaff = dim(R^"'"^) = 
3), we obtain 

Inserting this into (43), and observing that the cutoff rj eliminates the boundary 
terms at dKt^ , dKt^ , we obtain 



/ 4l*"V'a|' + pTapd'^d^r] dxdt = f TL0i^ + {d^r^Yp) da. 

•^^f[t2,«ii P -^a^fts,*!) P (47) 

Consider the second term on the left-hand side, which simplifies to pToarj" . This 
is supported on the regions t2 < t < t-ijl and 2t\ < t < ti, and we have the 
crude bound prj" = 0{l/t) on those regions. Thus by the energy bound (45), the 
contribution of this second term is 0{Eq). Now consider the right-hand side of 
(47). From (18) it is easy to show^"^ that T^^ = 0{Tlq) for all 13. Also, we have 
the pointwise bound 

\^ + (d''T^)p\<C{\e)<oo 
P 

on 9-K'(t2,ti), for some quantity C{\,e) depending on A and e (the value of C may 
vary from line to line). Thus we have 

/ ^\x"ipa? dxdt < CEo + C(A,e) / Tlq da. 

6) we thus have 
msup / ^ 

uniformly for all choices of A and e. We shall not take advantage of the ability to 
set e to be small, and just take e := 1. Then by the triangle inequality we have 

On the other hand, using the crude bound ^ < ||p and (45) we see that 



Applying flux decay (46) we thus have 

limsup / ■^\x°''4)a\'^ dxdt < CEq 



^l^il^l^ol^ dxdt < CEo 

, P 



"'^'^This is best seen using the null frame L, y, and L:=dt + dr\ indeed we have Tll = IV'iPi 
'^Lf = (V-i-V-y), and Tll = -^IV-yP- 



GEOMETRIC RENORMALIZATION OF WAVE MAPS 



25 



uniformly in A and ti, and thus 



limsup / ns l^'^l^ '^^'^^ - 

ti^O- J Kf^xH/2,2ti] 1*1 



In particular, for any A we see that for ti sufficiently close to (depending on A) 
we have 

l^^sf dxj<CEo. 



L 



'^[Ati/2,2ti] 

Telescoping this, we obtain that 

limsup ^ [ \\^s? dx^ < CEo/logX 

(note that one can replace the interval [— 1,T] by [to)?^] for any fixed to without 
affecting the limit superior). Letting A ^ oo we obtain (42). ■ 

Remark 3.2. The estimate works for all targets N, as it did not use the hypothesis 
that N had constant negative curvature. It is also completely covariant and did 
not rely the fixing of a gauge. 

Remark 3.3. Observe that in the above proof wc have in fact obtained stronger 
estimates than (42); for instance, we can improve the decay of jtps near the light 
cone, and obtain good estimates on relatively narrow slabs such as K[xt,t] instead 
of the tall slabs K^^i^x] in (42). For instance, it is not difficult to use the above 
arguments to obtain a sequence of times tk ^ 0~ , a sequence of scales A^ ^ oo, 
and Sk — > 0+, such that 

ml 



/ 

Jk, 



{p + ek\t\r 



dxdt < Ek- 



Morally speaking, this suggests that the rcscalcd original wave maps (j){t/tk, x/tk) 
are becoming locally self-similar as k oo. Unfortunately this is not by itself 
enough to deduce that these rescaled wave maps converge (weakly) to a non-trivial 
self-similar wave map, even after other rescalings, because the above estimates do 
not prevent the energy from dispersing into multiple points of concentration, or 
disappearing into the light cone; the basic problem is that control of the ipS compo- 
nent of the energy density docs not seem to directly control the other components 
without further structural control on <p, even near the light cone (the presence of 
angular components in the stress energy seems to prevent the Gronwall inequality 
approach in [7], [50], etc. from being effective, and naive attempts to exploit the 
negative curvature of the target N seem to require more boundedness control on <p 
than is currently available.) 
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